Demagnified GWs from Cosmological Double Neutron Stars and GW Foreground 

Cleaning Around 1Hz 
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Gravitational waves (GWs) from cosmological double neutron star binaries (NS+NS) can be 
significantly demagnified by strong gravitational lensing effect, and the proposed future missions such 
as BBO or DECIGO might miss some of the demagnified GW signals below a detection threshold. 
The undetectable binaries would form a GW foreground which might hamper detection of a very 
weak primordial GW signal. We discuss the outlook of this potential problem, using a simple model 
based on the singular-isothermal-sphere lens profile. Fortunately, it is expected that, for presumable 
merger rate of NS+NSs, the residual foreground would be below the detection limit Qaw,iim ~ 10~ 16 
realized with BBO/DECIGO by correlation analysis. 



I. INTRODUCTION 



A GW background from the early universe is one of 
the primary targets of observational cosmology. It would 
provide us with crucial information on physics at very 
high-energy scales. Among others, a primordial GW 
background generated during an inflationary period is a 
key objective [l[. Currently, there are two approaches to 
probe the inflation background. One is indirect observa- 
tion around / ~ 10 _18 Hz through B-mode polarization 
of cosmic microwave background Q. Another is direct 
GW detection around / ~ 0.1-lHz with the proposed 
space laser interferometers such as the Big Bang Ob- 
server (BBO) [|[ or the Deci-hertz Interferometer Gravi- 
tational Wave Observatory (DECIGO) 0,[1[. These two 
approaches at widely separated frequencies are compli- 
mentary and we might disclose fundamental properties 
of inflation by using them simultaneously (see e.g. 0). 
But, in both of them, we must cope with strong astro- 
physical contaminations to uncover the inflation back- 
ground. 

At present, the overall profile of astrophysical GW fore- 
grounds around / ~ 1Hz is unclear. Although white- 
dwarf binaries (important for the Laser Interferometer 
Space Antenna (LISA) Q) would not make a critical limit 
there @, we might have a strong foreground component 
whose quantitative properties are difficult to predict now. 
This is partly due to the complicated astrophysical pro- 
cesses involved. For example, to estimate the foreground 
made by supernovae of population III stars, we need, at 
least, the formation rate of these stars and their angular 
momentum distribution. But, unfortunately, these im- 
portant elements are poorly known at present Q . Mean- 
while, around IHz, we also have a foreground made by 
double neutron star binaries (NS+NSs) that are indi- 
vidually very simple system accurately predicted with 
the post-Newtonian inspiral waveforms [281 ] . In addition, 
NS+NSs are promising target for the first detection of 
GWs with ground-based detectors, and their coalescence 
rate has been extensively discussed with the observed 
abundance of NS+NSs in our galaxy [T2|. Considering 
our current understanding on the IHz band, in this pa- 



per, we study only the NS+NS foreground, neglecting 
other potential but highly uncertain astrophysical fore- 
grounds. 

NS+NSs can be regarded as solid GW sources around 
IHz, and cleaning of these binaries would be a critical ele- 
ment for the success of BBO/DECIGO to detect a weak 
primordial GW background. The basic approach for 
the cleaning is to identify individual binaries and sub- 
tract their chirping waveforms from the data of detectors 
As we see later, in order to fully use the designed 
sensitivity of detectors for a primordial background, the 
residual astrophysical foreground after the cleaning must 
be ~ I0 5 times smaller than their original strength (in 
terms of GW energy spectrum) . 

To discuss the prospect of this cleaning procedure, we 
first need to understand the detectability of GWs from 
individual NS+NSs. Since we cannot expect a large fluc- 
tuation (e.g. a factor of 2) for the intrinsic chirp mass dis- 
tribution of NS+NSs, the primary parameter that char- 
acterizes the signal-to-noise ratio of a binary would be 
its distance or equivalently its redshift z. Another im- 
portant parameter of a binary is its inclination angle, 
which can change the GW amplitude by a factor of 2y2 
(ratio between face-on and edge-on binaries). Therefore, 
the basic requirement for the cleaning is to make detec- 
tors that have enough sensitivities to detect an edge-on 
NS+NS at high redshift (e.g. z = 10) [nj. 

However, the situation becomes complicated due to the 
gravitational lensing effect that modulates observed am- 
plitudes of GWs during their propagation between the 
sources and detectors [131 ] . We might miss some demag- 
nified GW signals which are below a detection threshold, 
and their resultant residuals might be an obstacle for de- 
tecting a weak inflation background. Our principle aim in 
this paper is to provide a rough outlook on this potential 
problem caused by gravitational lensing. 

Gravitational lensing effects can be broadly divided 
into two categories: the weak lensing due to accumulated 
small distortions during wave propagation [l4| (see also 
[l5L [lt| for GWs from NS+NSs) and the strong lensing 
caused by specific massive objects with large distortions 
[HI]. We are interested in significantly (e. g. 50%) demag- 
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nified lensing fluctuations, but such probability is known 
to be negligible for weak lensing, even for a high-redshift 
source [17j . Therefore, we concentrate on the strong lens- 
ing effect that can generate significant demagnification 
more frequently. For the lens profile, wc use the singular 
isothermal sphere (SIS) model, which is fairly success- 
ful for studying various observational aspects of strong 
lensing effects [la ]. 

This paper is organized as follows. In Sec. II, we dis- 
cuss the merger rate of cosmological NS+NSs and their 
GW foreground. In Sec. Ill, the strong gravitational 
lensing effect is studied. The probability distribution 
function for the faint-end of demagnification is evaluated 
with the SIS model. Then we estimate the GW fore- 
ground made by undetectable NS+NSs. Section IV is 
devoted to discussions on this paper. 



II. COSMOLOGICAL NS+NS BINARIES 

In this section, we briefly discuss basic aspects of cos- 
mological NS+NSs and their GWs around 1Hz, follow- 
ing Cutler and Harms [io| . First we evaluate the to- 
tal merger rate Nt (T, total) of cosmological NS+NSs. 
Based on the observed NS+NSs in our galaxy and the 
number densities of various types of nearby galaxies, 
the comoving merger rate at present fiQ is estimated to 
be 10 _8 -10 _6 Mpc _3 yr _1 , roughly corresponding to the 
Advanced-LIGO detection rate 10 1=tl yr _1 0- To deal 
with merger events at cosmological distances, we put the 
comoving merger rate h(z) at redshift z as follows 



h(z) = s(z)fiQ. 



(1) 



Here, the nondimcnsional function s{z) represents the 
redshift dependence of the rate, and, in this paper, we 
consider the following two concrete models: (I) the fidu- 
cial evolutionary model sj(z) with a piecewise function 



1 + 2z (0 < z < 1) 
Sl ( z ) = { 3(5-z)/4 (l<z<5) 
(5 < z), 



and (II) a simple model 



sii(z) 



1 (0 < z < 10) 
(10 < z) 



(2) 



(3) 



with a flat merger rate up to z — 10. The first model s/ 
is the same as that used in Cutler and Harms [l(| ■ 

For estimating the total merger rate Nt observed to- 
day, we need to calculate the comoving volume element. 
The comoving distance to an object at redshift z is given 
by 



r(z) 



dz 



(4) 



where H{z) is the Hubble parameter at redshift z. In 
this paper, we fix the cosmological parameters at Hq = 



70km sec _1 Mpc _1 , Vt m = 0.3 and tt A = 0.7 for a flat 
universe. The function H(z) is written by 



(5) 



Since the comoving volume of a shell between z and 
z + dz is given by / LTrr(z) 2 (dr /dz)dz, the total merger 
rate is expressed as [To| 



dz 



r(z) h(z) 
(l + z)H(z) 



(6) 



with the redshift factor (1 + z) -1 due to the cosmologi- 
cal time dilution. Here, wc neglected a tiny increase of 
the actual number of merger signals due to the multiple 
lensed signals. For the fiducial model s(z) — sj(z), we 
numerically obtain 



N, 



1.1 x 10 5 



"o 



10- 7 Mpc- 3 yi- 1 



yr 



(7) 



For the flat-rate model s(z) — sjj(z), the prefactor in 
Eq.© becomes 5.9 x 10 4 . In Fig. 1, we show the relative 
redshift distribution of the merger event N(> z)/Nt- 
The numerator of this ratio is defined by 



N(> z) = Aire 



dz 



r(z) 2 h(z) 

(i + z)H( z y 



Next we evaluate the spectrum of the total GW en- 
ergy flGW.BT (B, binary; T, total) emitted by the chirp- 
ing NS+NSs before their subtractions. Here, following 
the standard convention, the energy density £Igw,bt is 
defined per the logarithmic frequency interval and nor- 
malized by the critical density of the universe. Its formal 
expression is given by [l8| 



n(z) 



... 8^/ 3 G 5 / 3 Af c 5/3 / 2/3 f°° J 

n GWtBT (f) ^jf2 j Q dz- + z)i/3H{z) . 

(9) 

In Eq.© the frequency / is measured at the detector 
frame, and M c is the intrinsic (not redshifted) chirp mass 
at the source frame. We keep using these two definitions 
throughout this paper. We also fix the chirp mass at 
M c = 1.22M©. The redshift dependence of the integral 
h(z)/[(l + z) 4 / 3 H(z)] in Eq.© can be decomposed into 
two factors: the event rate proportional to 

r{zfh{z)/[{l + z)H{z)] (10) 

as given in Eq. ^ , and the square of the individual signals 

h 2 oc l/Kz^l + z) 1 / 3 ] (11) 

(see discussions later in this section). For our fiducial 
model s(z) — si(z), we numerically obtain 



GW,BT 



(/) =4.0x10- 



n 



10- 7 Mpc- 3 yr" 1 



1Hz 



2/3 



(12) 
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FIG. 1: Fraction of the observed merger rate Nt and the 
foreground amplitude Qgw,bt by NS+NSs more distant than 
a given redshift. We plot numerical results for the fiducial 
merger rate s(z) = si(z) (see also [Tol]). 



while the prefactor becomes 2.1 x 10~ 12 for s(z) — su(z). 

The correlation analysis is a powerful method to de- 
tect weak stochastic GW signals [ijl]. In terms of the 
normalized energy density flawi we can i m principle, 
realize a factor of - SN^ifT^) 1 / 2 (SN, the signal-to- 
noise ratio for detection; T b s , observation period) better 
sensitivity, compared with the sensitivity ^GW,detector ~ 
10" 14 directly corresponding to the noise spectrum of 
BBO. With BBO/DECIGO and an observational period 
of T b s ~ 10yr, the limiting sensitivity of the correlation 
analysis to a stochastic GW background is SJc^iim ~ 
10~ 16 around the optimal frequency / ~ 0.3Hz where we 
are free from t he p otential foreground by double white- 
dwarf binaries 20]. This limiting sensitivity Qcwjim 1S 
a not a simple power-law function of GW frequency /, 
because of the shapes of the overlap reduction function 
(for definition see e.g. [Hj]) and the detector noise spec- 
trum. In order to fully exploit the potential specification 
of the proposed detectors and to pursue a weak primor- 
dial GW background down to their limiting sensitivity 
ClGW,lim ~ 10~ 16 , we need to subtract the NS+NS fore- 
ground and make the residual smaller by at least a factor 
of 



I / in 



5.6 x 10" 



BT 



»0 



10- 7 Mpc- 3 yr- 1 



(13) 



around the optimal frequency of BBO/DECIGO / ~ 
0.3Hz, by identifying NS+NSs up to their highest red- 
shift. Note again that the above target residual level 
(Q,GW,lim ~ 10~ 16 ) is much lower than the detector noise 
level (corresponding to VLaw,detector ~ 10~ 14 ). In Fig. 1, 
we show the fraction of the foreground &gw,bt made by 
binaries more distant than a given redshift z. 

Now we discuss the observation of chirping GWs from 
individual NS+NSs with BBO/DECIGO. Here, we study 
typical binaries using the unperturbed GW amplitudes. 
The lensed signals would be analyzed in Sec. III. In what 
follows, we neglect dependence on sky positions and po- 
larization angles of binaries, and we apply the averaged 



response of detectors with respect to these parameters. 
This is because multiple detectors with different orienta- 
tions would be used for BBO/DECIGO Hi], and we can 
expect a relatively weak dependence of signal-to-noise ra- 
tio on the direction and polarization angles, owing to 
an effective averaging effect (see e.g. |2l[ for averaging 
on the direction angles). This prescription simplifies our 
analysis (see also [l(|). But, with only one detector, more 
detailed studies would be required for sources with short 
signal durations (e.g. 1 week). 

We can evaluate the amplitudes of the two polarization 
modes of a binary with the quadrupole formula. In the 
frequency domain and in the principle polarization coor- 
dinate, their explicit forms for a circular orbit are given 



(h + ,h x ) = 



5 ^ 2 / 3 G 5 / 6 M c 5/6 /- 7 / 6 

(l + 2 ) 1 /6 c 3/2 r ( z ) 



(l+u 2 ,2u). (14) 



Here, wc defined the geometrical parameter u = cos / 
using the inclination angle /. A face-on (edge-on) binary 
has I = (/ = ±7r/2 respectively). Meanwhile, for a 
NS+NS at frequency /, the time Tqw before the merger 
is given by 



Tqw ~ 1 



/ 



0.2Hz 



-8/3 



(1 + *) 



-5/3 



yr. 



(15) 



and smaller than planned operation period of 
BBO/DECIGO. 

For signal analysis of each binary, we assume that, with 
some workable methods, the residual foreground would 
eventually become smaller than the detector noise af- 
ter subtractions of binary signals 10'. This justifies our 
treatment below in which we consider only the detector 
noises for estimating the signal-to-noise ratio of individ- 
ual binary. Note that this is a much weaker assump- 
tion compared with the previous requirement for realiz- 
ing the ultimate sensitivity for detecting the weak GW 
background with correlation analysis. We put the signal- 
to-noise ratio p(z, u) of unperturbed chirping GW from 
a NS+NS at a redshift z as follows 



p(z,u) = 



1/2 



BMc /6 F(u) 

(l + z)!/6 r ( z ) 



(16) 



with a function F(u) = 1 + 6u 2 



(see [13). Here, 



B is a constant determined by the noise spectrum of the 
detectors. Given the designed sensitivity of BBO, we fix 
the parameter B so that the signal-to-noise ratio becomes 
p(z = 5, u = 0) = 40 to match [l(|. For other redshifts, 
we have p(z = 10, 0) = 30, p(3, 0) = 52 and p(l, 0) = 110. 
We also have p(z, l)/p(z,0) = (F(l)/ '^(O)) 1 / 2 = 2^2 for 
dependence on the inclination angle, as commented in 
Sec. 1. For the fiducial merger model s — sj(z) with 



a realistic normalization n < 10 6 Mpc 



yr" 



Cutler 



and Harms [10( showed that the designed sensitivity of 
BBO [corresponding to p(z = 5, 0) = 40] would enable 
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us to detect all of the (unlensed) binaries and make their 
residual below the limiting sensitivity SJg^j^ ~ 10 -16 . 

For an evaluation later in Sec. IV, we also define the 
mean signal-to-noise ratios by 

_ _ 47re Jo°° dz Jq 1 du ggffg) p(z, u) 2 

We numerically evaluated this expression and obtained 
p = 187 for s = sj and p = 161 for s = sjj. 



III. DEMAGNIFIED GW SIGNALS 

A. Lensing Probability 

Strong gravitational lensing produces multiple signals 
for an intrinsically single event 13]. In this subsection, 
we discuss the faint end of GW signals caused by strong 
lensing. For the density profile of the lens objects, we use 
the SIS that captures lens structures well at the scales rel- 
evant for the strong lensing events [l3j |. Here, we do not 
include detailed effects such as the external tidal shear 
field or the ellipticity of lenses, which cause minor cor- 
rections for the strong lensing probability itself (see e.g. 

[23l ] . and discussions on the quadruple lenses therein). 

2 

The explicit form of the SIS density profile 37r v GR 'i 
(R: the radial coordinate) is characterized by the one- 
dimensional velocity dispersion v. We further introduce 
the unperturbed angular diameter distances D OSl D i and 
Di s between the observer-source, observer-lens and lens- 
source respectively. For example, the distance D a i is 
given by 



Dal = 



r(zi) 



1 



zi 



(18) 



with the comoving distance r(z) defined in Eq.(j4|. In 
order to study lens mapping, we define the characteristic 
angle 



9 d = 



Attv 2 D 



1.1 



c 2 D os 



(19) 



and use the two-dimensional coordinate y for the source 
plane normalized with the length unit 6dD os . We also 
introduce the coordinate x for the lens plane normalized 
with the length BaPoi- The origin of these two coordi- 
nates are fixed at the direction of the lens center from 
the observer. Then the lens equation for the SIS model 
is given by [l3[ 



y = x — xj \x\ 



(20) 



Because of the apparent symmetry around the lens cen- 
ter, the lens mapping between x and y is essentially one- 
dimensional correspondence on lines passing their origins. 
The lens equation ([2U|) has two solutions 



x± = y ± 1 



(21) 



for a small impact parameter at < y < 1, and one 
solution 



y - 



1 



(22) 



for a large impact parameter at y > 1. 

The amplification factor A of gravitational lensing is 
evaluated with the Jacobian of the mapping and given 
by 



A = 



x dx 
y dy 



(23) 



For the two solutions x± , the amplification (f2"3")l becomes 



A, 



y/T/y + T, A_ = y/l/y-1. 



(24) 



Comparing these two signals, the second one corresponds 
to the inner, fainter signal with a later arrival time. Note 
that the amplification A is defined by the square root 
of the magnification p which is often used in the liter- 
ature on gravitational lensing of electro-magnetic waves 
[l3l | . Our definition here reflects the important fact that, 
for detecting GW from a binary, we directly observe the 
waveform itself (oc h), not its energy (oc h 2 ). 

In contrast to the identity A+ > 1, the fainter coun- 
terpart A_ approaches zero in the limit y — > 1 (from 
below). Since we are interested in weak GW signals that 
might be undetectable, we mainly study the fainter one 
A- rather than the brighter one A+. The signal-to-noise 
ratio of a demagnified signal is expressed as A-p(z, u). 

At y ~ 1, the time delay between two signals is given 

by m 



At 



00 2 m 4 d o1 d 1s 

327r y-J — + 



cD os 



8.6 x 10 b 



210km/sec 



D ol D ls D-}(l + Zl ) 
lGpc 



sec. 



For a typical velocity dispersion v of a lens object, we 
have the relation fAt > 10 6 3> 1 around the optimal fre- 
quency of BBO/DECIGO / - 0.3Hz. Therefore, in the 
present situation, the geometrical optics approximation 
would work well [2J] , and the strongly- lensed two signals 
would be observed as two distinct chirping waveforms. 
But it might be interesting to analyze potential wave ef- 
fects (e.g. in relation to substructure of lenses) (24j. For 
the SIS model, the separation angle A9 between the two 
images is given by 



AO = 28 d - 2.4" ( — » 1 



D os J \210km/sec 



(25) 



With the expression for the characteristic angle 84 de- 
fined in Eq. (fT9| . the probability of the strong lensing 
event for a source at redshift z s is given by 



Pd(Zs) 



dzi 



dv 



dN 167r 3 C -V(l + zi) 2 Df s D 2 ol 



dv 



H{zi)Dl s 



(26) 
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where is the distribution function for velocity disper- 
sions of lens objects. In this paper, we use the following 
model 



dN 



<t>*P 



— exp 



(27) 



characterized by parameters 4>* = 4.1 x 10 3 /i 3 Mpc 3 , 
w* = 88.8km sec" 1 , a = 6.5 and /3 = 1.93 [25]. We 
neglect the cosmological evolution of this function. This 
is partially supported by actual lensing observation up 
to z s ~ 1. Our main results [e.g. Eqs. ([37f and (|H|) ] 
change only slightly, for an evolutionary model with the 
following set of parameters analyzed in [25] 



(j>*{z{) = (j)*(l + Zi) 



-0.229 



, (T*(zi) = CT*(1 + Z() 



-0.01 



(28) 



The product v^dN/dv is relevant for the lensing proba- 
bility and also for the distribution of time delays. It be- 
comes maximum at v — 2.4er» = 210km/sec and steeply 
declines at large v, as expected from the fraction 



/ 4 ~ dvv 4 dN/dv 
J °° dvvHN/dv 



10" 



(29) 



Even for a high-redshift source z s ~ 5, the redshift inte- 
gral (|26p has dominant contribution around zi ~ 1. In 
Fig. 2, we present the lensing probability Pd(z), as a 
function of the source redshift z. 

On the source plane, the demagnification A- (< 1) cor- 
responds to the normalized radial coordinate as 



1 



y 



1 



(30) 



Therefore, the probability that a source at a redshift z 
has a demagnification in the range [A—, A- + gL4_], is 
given by 



dP 

HaZ 



dA- = P d {z) 



4A- 



(A 2 _ + 1) 



-dA^ 



(31) 



B. Undetectable Lensed Signals 

In this subsection, we study the subtraction problem 
for lensed GWs emitted by cosmological NS+NSs. First 
we evaluate the event rate for strongly-lensed GWs. For 
binaries at redshift z, the fraction Rq(z) of the strongly- 
lensed GWs coincides with the probability Pd(z) given in 
Eq.flHJ) as 



Ro(z) = P d (z). (32) 
Then the total merger rate of the lensed signals is written 



Nr. 



4-7T 



dz 



r(z) 2 R (z)fi(z) 
(l + z)H(z) ■ 



(33) 
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FIG. 2: Probability of having multiple lensed images as a 
function of the source redshift z. For the velocity distribu- 
tion function dN/dv, we use the no-evolution model given by 
Eq.(H3. 



For the fiducial model s = Si{z), we numerically evalu- 
ated this expression and obtained the averaged lensing 
probability 



P = N L /N T = 3.7 x 10" 



or equivalcntly 



Nr = 40 



»0 



10- 7 Mpc- 3 yr~ 1 



yr 



(34) 



(35) 



For the flat-rate model s(z) = sn(z), we have the ratio 
N L /N T = 8.6 x 10~ 4 . 

Next we discuss a demagnified GW signal whose signal- 
to-noise ratio A_p(z, u) is less than a detection threshold 
Pth- To begin with, we estimate the relevant value pth, 
taking into account the observational situation of the two 
lensed signals. The typical time delay At between the 
earlier bright signal with A+ and the later faint one with 
A_ is much less than the planned observational period 
T obs - lOyr of BBO/DECIGO. Thus we mainly consider 
a scenario to search for the demagnified signals under the 
detection of the bright ones [29j. Except for the coales- 
cence time t c , the fitting parameters of the normalized 
waveforms of the two lensed signals can be regarded as 
almost identical, and these parameters would be gener- 
ally well determined with the bright first image. For 
example, with BBO, the typical size of the localization 
error-ellipsoids in the sky is ~ 10" x 10" for a NS+NS at 
z ~ 5 10]. On the other hand, the characteristic image 
separation AO in Eq. pS)) is several arcseconds. There- 
fore, directions of two lensed images would be fitted by 
the same parameter for a NS+NS at z ~ 5. But, for a 
NS+NS at low redshift, the second image could get out 
of the error-ellipsoid of the brighter image. 

The similarities of fitting parameters would signifi- 
cantly ease the detection of the demagnified one with 
the matched filtering method, compared with its sin- 
gle detection for which we need 10 30 -10 36 templates to 
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make a full coherent signal integration. Even with the 
estimated computational power available at the era of 
BBO/DECIGO, we cannot make so many templates, and 
we need to use a suboptimal detection method which re- 
quires a higher detection threshold (p cr > 20 see (Tpj ) 
than the (optimal) coherent integration. Therefore, the 
information of the brighter signal considerably decreases 
the detection threshold p t h for a lensed second signal, 
compared with its single detection. For the typical num- 
ber of templates n tc ~ 10 7 -10 8 estimated from the bin- 
ning of the coalescence time t c of the faint image, the 
required condition for the threshold pth with respect to 
a false alarm rate f r is given as 



fr ~ N tc eri(p th /V2), 



(36) 



and we have the solution p t h ~ 6.5 for f r ~ 0.01. Includ- 
ing a safety factor (e.g. potential resampling of other 
fitting parameters around the bright signal), we take 
Pth = 10 as a standard value for detecting the faint sec- 
ond image. For a reference we also use a pessimistic value 
pth ~ 20 in the analyses below [Io| . The results for this 
higher value would also provide us with an insight about 
the stand-alone analysis for a demagnified signal without 
using the information of the associated brighter one. 

Here, we briefly comment on the identification of lensed 
pairs. The order-of-magnitude estimation for the re- 
quired numbers of binning for the redshifted chirp mass 
M c (l + z) is ~ 10 8 , which is much larger than the relevant 
numbers of NS+NSs ~ 10 5 [10]. In addition, given the 
good localization in the sky, directions of binaries would 
also become useful information to identify a lensed pair. 
Therefore it is unlikely to misidentify GWs from two dif- 
ferent NS+NSs as a lensed pair. Although our aim in 
this paper is not to discuss further scientific possibilities 
with identified lensed pairs, such studies would be also 
interesting [l5|, 

EH- 

Now we statistically study the demagnified GW signals 
that are below the detection threshold with A_p(z,u) < 
Pth- The fraction of the number of undetectable signals 
at a redshift z is given as follows; 

SfdA-fidufceQHH-A-pjz,*)) 
Ri(z) ee -5— . (37) 

Jo du 

Here Q(x) is the step function. In Fig. 3, we plot the 
function R\(z) for two choices: p t h = 10 and 20. At 
the small amplification regime A- -C 1, the probability 
of having an amplification less than a given value A_ is 
proportional to A 2 _ [see Eq.(l3"Tj)], and we have an asymp- 
totic scaling relation R\ oc pf h . But we should be careful 
to note that this scaling relation depends on the details 
of the lensing profile and is not universal. The event 
rate Nn (U, undetectable) of the undetectable signals is 
evaluated as in Eq. (f33|) and we have 



Nr. 



An 



dz 



r(z) 2 i?i(z)n(z) 



For the fiducial model s(z) — sj(z) and the threshold 
Pth = 10, we numerically obtained the result 



Nr, = 1.1 



7f() 



10- 7 Mpc- 3 yr- 1 



yr 



or equivalently 



Nu/N T = 1.1 x 10" 



(39) 



(40) 



With the higher threshold pth = 20, the ratio becomes 
Nu/N T = 3.9 x 10~ 5 . 

We move to evaluate the energy spectrum of the GW 
foreground made by undetectable binary signals with 
A-p(z,u) < pth- Here we introduce the factor R2{z) as 
the fraction of the GW energy due to the undetectable 
signals at redshift z relative to the total binaries at the 
redshift. Summing up the contribution of undetectable 
ones, the factor R2(z) is formally given by 



R2(z) = 



(38) 



fdA. ^ du^-e(pth - A_p(z,u)){A_p(z,u)} 2 

Jp 1 dup(z, u) 2 

(41) 

In Fig. 3 we provide the function R2(z) for two thresholds 
Pth = 10 and 20. In the present case, we can derive the 
asymptotic behavior R 2 oc p\ h for the SIS profile, after a 
simple consideration on the factors in the numerator of 
Eq.gl]). 

With the fraction R 2 (z), the energy spectrum flaw Fin 
of the undetectable demagnified ones is expressed as [30| 

_ 8^M 5 J 3 p/z p n{z)R 2 (z) 
nGW ' BU= 9B2 J (l + z)^H(zY (42) 

With the fiducial model s(z) = si(z), we numerically 
obtained SIgw,bu /&GW,BT — 1.5 x 10~ 8 for p t h = 10 
and &gw,bu /^gw.bt = 2.1 x 10~ 7 for p th = 20. With 
s(z) = Sjj(z), we have 9.8 x 10 -8 (p th = 10) and 
1.2 x 10~ 6 (pth = 20). Compared with the ratio in flS]). 
the undetectable GW foreground would be comfortably 
smaller than the limiting sensitivity £lGW,Um ~ 10~ 16 for 
the realistic normalization ho < 10 _6 Mpc _3 yr _1 . 



IV. DISCUSSIONS 

So far, we have studied the demagnified GW signals 
with the SIS density profile. Although this lens model 
is quite simple, it has reproduced observational results 
(with electro-magnetic waves) of gravitational lensing 
fairly well [13| . But one might think that this model 
has been discussed at the relatively high- amplification 
regime where observations can be performed more easily. 
In this respect, it might be reasonable to wonder whether 
the simple SIS model is also a useful tool to study the 
strong lensing effects at the low-amplification regime as 
discussed in this paper. One simple example highlighting 
the situation is about the central third images of strong 
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to this example, it should be mentioned that, compared 
with typical electro-magnetic wave observation, GWs are 
generally measured at much lower frequencies and could 
be more susceptible to the wave effects (e.g. by sub- 
structures) that could complicate the signal analysis [24[ • 
In relation to the foreground cleaning which is essential 
for directly detecting a GW background from inflation, 
further studies on the faint end of lensing amplification 
would be worthwhile beyond our simple treatment. 

On a final note, we attempt a more robust approach 
to evaluate the strength of the foreground composed by 
residual GW signals from undetected binaries, without 
using details of the probability distribution (13~T1) for the 
demagnification A—. Here, we simply assume that, if 
GW is strongly lensed, then the total signal-to-noise ratio 
of its unsubtracted component would be less than the 
detection threshold pth- We define the residual Qgw,bs 
(S; simplified) as the summation of these unsubtracted 



ones. Since we have the relations Q 



GW.BT 



oc N' 



TP' 



and 



Qgw,bs oc N L p th 2 = N T Pp th 2 (see Eq.flM])), we obtain 



GW,BS 



P 



GW.BT 



Pth 
P 



(43) 



We regard the amplitude flGW,BS as a conservative up- 
per bound for the potential foreground made by the un- 
detected signals. With the numerical results presented 
so far, we can evaluate this ratio and obtain 



FIG. 3: Relative contributions of the undetectable GW sig- 
nals from NS+NSs at redshift z. The BBO noise curve is as- 
sumed with two choices of the detection threshold pth (dashed 
curve: p t h = 20 and solid curve: pth = 10). Left panel: 
the number fraction Ri(z) of the undetectable signals. Right 
panel: the fraction Ra(z) of the GW energy density due to 
the undetected NS+NSs. 



lensing [26|, |27j ■ A lensing galaxy could have a core struc- 
ture around its center, instead of the power-law profile 
cx r -2 of the SIS model, and such a density profile can 
generate a very faint third image around the direction 
toward the center of the lens. However, observational 
analysis of the third images is known to be q uite difficult, 
partly due to the faintness of the signals [3l| . In addition 



GW,BS 



= 1.1 x 10" 



GW.BT 



(44) 



for the fiducial model s(z) = si(z) and 3.3 x 
10 _6 (p t / l /10) 2 for s(z) = sn(z). Therefore, unless the 
normalization ho is relatively high no > 10 _6 Mpc _3 yr~ 1 , 
the upper limit Qgw,bs is smaller than the limiting sen- 
sitivity flcwMm [see Eq. ([T5j) ]. and we expect that the 
residual foreground caused by demagnified signals would 
not be a fundamental problem to detect an inflation back- 
ground of fl G w > 10 -16 with BBO and DECIGO. 
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